Dissipation through spin Coulomb drag in electronic spin dynamics 
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Spin Coulomb drag (SCD) constitutes an intrinsic source of dissipation for spin currents in metals and semi- 
conductors. We discuss the power loss due to SCD in potential spintronics devices and analyze in detail the 
associated damping of collective spin-density excitations. It is found that SCD contributes substantially to the 
linewidth of intersubband spin plasmons in parabolic quantum wells, which suggests the possibility of a purely 
optical quantitative measurement of the SCD effect by means of inelastic light scattering. 
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Spintronics applications are receiving increasing attention 
in the hope of revolutionizing traditional technology by a pow- 
erful exploitation of the spin - as well as the charge - de- 
grees of freedom. An intense research effort is underway to 
improve our understanding of spin dynamics, especially re- 
lated to nanocircuits and their components, such as quantum 
wells and wires. In this context the theoryof spin Coulomb 
drag (SCD) was recently developed Q B B [3 This the- 
ory analyzes the role of Coulomb interactions between dif- 
ferent spin populations in spin-polarized transport. Coulomb 
interactions transfer momentum between different spin popu- 
lations, so that the total momentum of each spin population 
is not preserved. This provides an intrinsic source of fric- 
tion for spin currents, a measure of which is given by the 
spin-transresistivity [1]. SCD is generally small in metals, 
due to a typical Fermi temperature of the order of 10 5 K, 
but can become substantial in semiconductors, where the spin- 
transresistivity can be larger than the Drude resistivity ylldl- 
As the quest for defect-free materials with longer and longer 
spin-decoherence times is continuing, spurred by practical re- 
quirements in spintronics as well as in quantum computation 
devices, the SCD is bound to become one of the most serious 
issues in spin polarized transport, since, due to its intrinsic 
nature, it cannot be avoided even in the purest material. In 
factjthe recent experimental observation of SCD by Weber et 
al. la] shows that the effect dominates spin diffusion currents 
over a broad range of parameters, in agreement with theoreti- 
cal predictions |2j,|3|,|5|]- 

In this letter we discuss a critical issue for potential spin- 
tronics devices, namely the power loss in spin transport and 
dynamics due to SCD. We shall analyze in detail its effect on 
optical spin excitations, and propose an experiment to mea- 
sure the intrinsic SCD linewidth enhancement of spin plas- 
mons in parabolic semiconductor quantum wells. While up to 
now SCD has been considered only in relation to spin trans- 
port, the proposed experiment would provide an alternative 
way of measuring this subtle effect, and thus establish un- 
equivocally the influence of SCD on optical excitations. 

Let us consider a system composed of spin-up and spin- 
down electron populations, as for example the electrons in the 
conduction band of a doped semiconductor structure. We are 
assuming spin-flip times long enough so that spin populations 



are well defined on the relevant time scales. This assumption 
- at the very core of spintronics - has been proved reasonable, 
with experimentally measured spin-decoherence times of the 
order of microseconds |7]. Previous papers on SCD have 
mainly analyzed the dependence of the spin-transresistivity 
over temperature this letter will focus on its fre- 

quency dependence 1 1 ], which is important both for AC spin- 
tronics applications and spin-resolved optical experiments. 

In the linear response regime and for weak Coulomb cou- 
pling one can write a phenomenological equation of motion 
for the spin a population 1 1]. The SCD force is defined as the 
Coulomb force (per unit volume) exerted by spin cr(= — a) 
electrons, moving with velocity v ff , on spin a electrons, mov- 
ing with velocity v CT : 



F CTff (r; uj) = -7(w)m- 



[v CT (r) - v ff (r)] 



(1) 



where the number density, n a , of er-spin electrons of effec- 
tive mass to, and the total density, n = + rij., are those 
of a homogeneous reference system. The drag coefficient 7 
appearing in Eq. Q is directly proportional to the real part of 
the spin-transresistivity /0|x(lj]: 



3fym( w ' T ; n T' r H) 



(2) 



where T is the electronic temperature. Sip-fj. has a nega- 
tive value and can be defined by the relation Ef|j _ = 
— ejip^i, with jcr the number current density of the a spin 
population, Ej the effective electric field which couples to 
the f-spin population and includes the gradient of the local 
chemical potential, and e the absolute value of the electronic 
charge. 

As noted above, SCD provides an intrinsic decay mech- 
anism for spin-polarized currents, and is thus a source for 
power loss in a spintronics circuit or device. From the gen- 
eral definition of power and using Eq. <□, the SCD power 
loss density per unit time for the er-spin population is given by 



P CT (r;o;,n T ,n;) = F^r) • v CT (r) 



|ja(r)| 2 - j ff (r) • j CT (r) 

TlfT 



(3) 



(4) 
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FIG. 1: Spin-transresistivity vs rescaled frequency Tllj/Ef 

for n = 10^ cm -3 , a; = 16, 17, 18 as indicated, and GaAs parame- 
ters (m = 0.067m e , e = 12). Inset: |5Rpn | in mf2 cm vs hui in eV 
for the same densities. Dashed line: high-frequency limit, Eq. |6|. 



Notice that P a can change sign depending on the relative 
strength and direction of the spin-resolved current densities, 
a positive sign implying that the a spin population is being 
dragged along by the faster a spin population. In particular, 
for a system with spin populations drifting at the same aver- 
age velocity, P a (v; to) = 0. The total power loss per unit time 
in a system with a slowly varying density can be calculated as 



d 3 r[P a {r;uj,n 1 {r),n l (r))} 



(5) 



Fig. [^shows the transresistivity $tp-\i (oj; , ni) as a func- 
tion of frequency, calculated for GaAs at T = 0, using a gen- 
eralized random phase approximation Q|. We see that $tp-\i 
has a maximum when Ep cr (n C r(z)) is of order fiu (Ep a is the 
cr-spin Fermi energy). This maximum roughly scales as [3] 
(ha*/e 2 )/n s « 140 pft cm • em e /(mn s ) with s<l: it is then 
reasonable to expect a sizable damping effect due to SCD. We 
notice also that for very low densities, i.e. Ep <C fid), 
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(6) 



independent of the carrier density (see Fig. IT]inset) lll4ll . 

Due to problems with electrical injection |8] and the neces- 
sity of driving spin dynamics on sub-picosecond time-scales 
tejl, large attention has been focused on optical spin injection 
H and optically controlled spin dynamics 1 10]; in the fol- 
lowing, we will explore how the SCD affects the lifetime and 
dynamics of spin-dependent optical excitations. 

The excitation spectrum of a system can be calculated in 
principle exactly with time-dependent density-functional the- 
ory (TDDFT) O- In TDDFT, the properties of an interact- 
ing time-dependent many-body system are described through 
a non-interacting time-dependent system (the so-called Kohn- 
Sham system) characterized by an exchange-correlation (xc) 
potential. The xc potential is a functional of the current I12B . 
and needs to be approximated in practice. An approximation 



which is cast in the language of hydrodynamics, including dis- 
sipation effects, was proposed in Ref. 11311 : nonadiabatic xc 
effects manifest themselves as viscoelastic stresses in the elec- 
tron liquid, which are proportional to the velocity gradient. 
The corresponding expression for spin-dependent systems is 
discussed in 1 14], the main difference being the appearance of 
a term - in addition to the viscoelastic tensor - representing 
the damping of the spin currents due to the SCD effect. 
Our derivation of the excitation energies for a spin-depen- 
systei 

m 

response equation, 



dent system closely follows the spin-independent case, see 
Ref. I ll 51 1 for details. Starting point is the TDDFT current 



j<x(r,w) 



(7) 



Here, x <x( r i r 'i w ) i s me Kohn-Sham current-current response 
tensor, which is diagonal in the spin-channel. The effective 
vector potential is defined as a CT = a° xt + a^ 1 + a* , where 
a™* is an external perturbation, and the Hartree and xc vector 
potentials are given by 



M 2 



(8) 
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iujn a (r) 



W )/ J ££,(r,r / ) 



— n <7 (r)n rf (r)p n {uj;n a (r)n s (r)) 



7i cr (r)n (T '(r) 
where v : k are Cartesian indices. In Eq. (|9j 



(9) 



drier drier' 



"ot,i(r) 
(10) 

is the frequency-independent xc kernel associated with the 
adiabatic local-density approximation (ALDA), where e£ c is 
the xc energy density of a homogeneous electron gas, and no a 
the ground-state spin density of the system. The other terms in 
Eq. (|9jl represent non-adiabatic xc contributions, which bring 
in the dissipation. In the second term, er™ aa , is the spin- 
resolved viscoelastic stress tensor of the electron liquid 
The key quantity in the last term of Eq. (|9jl is /9fj_. 

We now consider a specific excitation pa — > qa between 
the Kohn-Sham levels ip pcr and ifjqa, and assume the ground 
state to be spin unpolarized. To derive the TDDFT correction 
to the bare Kohn-Sham excitation energy Tiuj^a, we apply the 
so-called small-matrix approximation 1 15l fl7ll . The result is, 
to lowest order in the non-adiabatic corrections, 



fc2 2 



*2 2 
n 0J 



'pqcr T 2SW : 



ipqt 



\( QH + ALDA I CH+ALDA\ 

" IA°erer 31 D Sa ) 

f (s s ™±srj>)}, 



(11) 
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where the +/— sign refers to charge- or spin-density excita- 
tions (CDE/SDE) respectively. The terms in square brackets, 
S^+ ALDA , S™, and S*°?, are the dynamical many-body cor- 
rections to the bare transition energy fno vqa between the single 
particle levels pa and qa. The Hartree+ALDA shift is given 
by 



QH+A 



d 3 r I dV V>p<x( r )Vv ( r )VV' ( r ')</V ( r 



r — r 



_ , MLDA/ /N 
/ I ~ ./xc,<t<t' \ L 1 L / 



(12) 



which causes no dissipation, /* C L ° A / being frequency indepen- 
dent and real. The viscoelastic shift is given by 



= 



luJ 



E 



n CT (r) 



(13) 

where &%f£ a , is the xc stress tensor with the exact 

current \ a . y replaced by j P9CT (r) = (ip pa \3a\ipq<7), with j ff the 
paramagnetic particle current density operator. Eq. Jl 3I > can 
be viewed as the average rate of energy dissipation per unit 
time in a viscous fluid 1 15, 16], where cr^'^L is the viscoelas- 
tic stress tensor of the fluid, and V^'pga^/io-] me velocity 
gradient. In contrast to the familiar expression from classical 
fluid dynamics [16], 5 VE has both real and imaginary part. 
The SCD shift is a central result of this letter: 



riSCD i crSCE 



ie 2 to 



"pqcr 

n s (r) 
n a (r) 



rf 3 rpu( w ; n T( r ). n i( r )) 



iW(r)| 2 TWO) -w(r) 



(14) 



As we will show in an example below, under certain circum- 
stances this new contribution to the broadening of an excita- 
tion can actually dominate the damping process. 

By comparison with Eqs. and l|5}, we immediately rec- 
ognize the structure of the power loss typical of the Coulomb 
drag force 1 18]. Like the viscoelastic term (II 31 . the SCD term 
(114-i contains both a real and an imaginary part. Notice that, if 
the external driving force couples in a different way to the two 
spin components, such that the average spin velocities are dif- 
ferent, the SCD term contributes to the charge channel too. In 
this particular case the two spin-populations may be consid- 
ered distinguishable, characterized by a spin-dependent fre- 
quency both in the charge and in the spin channel. This 
implies that the Coulomb drag force exerted by one popula- 
tion onto the other can be regarded as an external force. 

This concept can be clarified by considering the charge and 
spin plasmons in a quantum well 11191 l20t 12111 . The inset to 
Fig. |2] illustrates the two types of density oscillations for a 
parabolic well, in which the rtf and components move back 
and forth in phase (CDE) or with opposite phase (SDE). In the 
case of the SDE, the average net momentum transferred by 
Coulomb interactions from the a to the er-spin population will 
be directed opposite to the er-spin direction of motion, so that 
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FIG. 2: Upper panel :Spin-plasmon linewidth for a parabolic 

quantum well versus curvature parameter A, for N s — 10 , 10 
and 10 12 cm" 2 and GaAs parameters. The inset illustrates the col- 
lective motion of the two spin populations (CDE: in phase, SDE: out 
of phase). Lower panel: rescaled linewidth Fg^/hcj versus A for 
the same system and parameters as upper panel. 



the SCD effect damps the motion of both spin populations. 
For the charge plasmon the effect can become more subtle: 
since the average spin velocities are in the same direction, the 
net result of Coulomb interactions between the two spin pop- 
ulations will be to transfer momentum from the "hotter" to the 
"colder" population, until equilibrium is reached. In this case 
the SCD effect would not damp the motion of both spin pop- 
ulations, but pump momentum from the faster to the slower. 

We now proceed to estimate the size of the SCD effect for 
optical excitations in a parabolic quantum well. According to 
the Harmonic Potential Theorem 1 22], the intrinsic linewidth 
of a CDE in a parabolic confining potential is strictly zero. 
The TDDFT linear response equation satisfies this require- 
ment: CDE's in a parabolic well have a uniform velocity pro- 
file, so that the viscoelastic stress tensor vanishes. Likewise, 
in expression Jl 31 for S™,, V K \j P qa,v/n a ] is very small. The 
viscoelastic contributions to SDE's are thus a higher-order 
correction compared to the SCD contributions, which give 
the dominant correction to the excitation frequency beyond 
ALDA. The intrinsic SDE linewidth for a parabolic quantum 
well therefore becomes Tsde 



r|gg, where 



^SCD 
SDE 



SCD 



e 2 N s u) 

^pqa 



rrSCD] 
(7(7 \ 

dz dip n {ui;n ] (z),n l (z)) 



n (7(z) ,. , V 2 
n a {z) l3pqa(Z) 



Tjpgo'(^) ' jpqa{ z ) 



, (15) 



with N s the two-dimensional electronic sheet density. 

Numerical results for for a GaAs-based quantum well 
are shown in Fig. [2] We assume only the first subband to be 
occupied, i.e., n a (z) — N s \il>i a (z)\ 2 , and approximate the 
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Kohn-Sham orbitals ip qtPa (z) entering Eq. dl5l by the first 
two eigenstates of a harmonic oscillator with external poten- 
tial Tv 2 z 2 /2m\ A . Furthermore, to lowest order in the non- 
adiabatic corrections oj a can be replaced with u) pqcr . For this 
system the parameters which govern the linewidth of the SDE 
mode are N s and the quantum well curvature parameter A. 
The latter determines both the excitation frequency and the 
characteristic width of the ground-state density distribution. 
The results in Fig. |2]show that r|§§ can be nonnegligible (a 
large fraction of me V) for experimentally reasonable parame- 
ters, and T%^/fiuj can be of the order of few percents for a 
large range of curvature parameters and carrier densities. 

For a specific N s , the linewidth exhibits a well defined max- 
imum as a function of A. The position of this maximum is 
determined by the competition of two distinct effects: (i) The 
low-density saturation value of p-|-j_ increases with A [i.e. de- 
creases with u, see Eq. 0]; (ii) The average particle velocity 
decreases with A (i.e. decreases with the parabolic curvature). 
The two effects give opposite contributions to the dissipation 
[see Eq. (0], and the maximum occurs when the second effect 
takes over. Due to the density dependence of p-pj_ (see Fig.Q, a 
substantial contribution to the integrand in Eq. (1151 can come 
from the lateral regions of the quantum well, where the parti- 
cle density is low. This is in contrast to the VE contribution, 
which tends to be dominated by the high-density regions. 

The above example shows that, even when other forms of 
damping, such as disorder and phonons, are drastically re- 
duced by careful selection of the system characteristics, the 
dissipation induced by SCD cannot be avoided, due to its in- 
trinsic nature. 

Eq. H5\ suggests an experimental way to extract the im- 
pact of SCD on spin dynamics by the optical measure of 
the linewidth of both charge- and spin-plasmons in the same 
parabolic quantum well. Such measurements can be carried 
out using inelastic light scattering |23[. Under the reasonable 
assumption that (i) extrinsic (ext) damping (non-magnetic im- 
purities, phonons) affect the CDE and SDE in the same way, 
and (ii) the viscoelastic term can be disregarded due to the 
parabolic system geometry, we have 

1 SDE — f CDE ~ (i SDE + 1 SDE) ~ i 1 CDE) ^ SDE ) ( ife ) 

i.e., the SCD contribution to the spin-plasmon linewidth is 
given to a very good approximation by the difference of the 
SDE and the CDE linewidths. This provides a valuable op- 
portunity for comparison with microscopic models for the 
transresistivity via Eq. J15I . using the appropriate Kohn-Sham 
single-particle orbitals of the system. 

In conclusion, we have presented a discussion of the power 
loss in a device due to dissipation of spin-dependent currents 
induced by SCD forces. We have suggested a new, purely op- 
tical method to measure the SCD effect in spin-density excita- 
tions in parabolic quantum wells. In the u> — > limit, a partic- 
ularly interesting application of our formalism would be to de- 
scribe the SCD intrinsic dissipation in spin-dependent trans- 
port through single molecular junctions ll24ll . As the broad 



effort in spintronics, quantum computation and transport in 
micro- and mesoscopic systems continues, we expect a grow- 
ing impact of the SCD effect in future applications. 
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